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ABSTRACT 

According to one version of the recently proposed "manifold" theory that explains the origin 
of spirals and rings in relation to chaotic orbits, galaxies with stronger bars should have a higher 
spiral arms pitch angle when compared to galaxies with weaker bars. A sub-sample of barred- 
spiral galaxies in the Ohio State University Bright Galaxy Survey, was used to analyze the spiral 
arms pitch angle. These were compared with bar strengths taken from the literature. It was 
found that the galaxies in which the spiral arms maintain a logarithmic shape for more than 70° 
seem to corroborate the predicted trend. 

Subject headings: galaxies: kinematics and dynamics — galaxies: spiral — galaxies: structure — 
galaxies: kinematics and dynamics 



1. INTRODUCTION 

Spiral arms in barred galaxies have been ex- 
plained in the past as density waves (e.g., Kor- 
chagin & Marochnik 1975) or spiral waves that 
result from the crowding of gas orbits (Huntley et 
al. 1978). Kaufmann & Contopoulos (1996) in- 
voked for the first time the need for chaotic orbits 
as building blocks of spirals near the end of the 
bar. In the Kaufmann & Contopoulos (1996) mod- 
els, regular orbits dominate the main structure of 
the bar and the outermost portions of spiral arms. 
The inner portions of spiral arms are supported by 
chaotic orbits. Recently it has been proposed that 
chaotic motion can support the spirals in barred- 
spiral systems. The new theory proposes that un- 
stable Lagrangian points (Li or L2) near the end 
of the bar are the sites where chaotic orbits are 
guided by invariant "manifolds", and are the ori- 
gin of spirals and (inner and outer) rings (Voglis & 
Stavropoulos 2006a; Patsis 2006; Romero-Gomez 
et al. 2006; Voglis et al. 2006b, c; Romero-Gomez 
et al. 2007; Tsoutsis et al. 2008, 2009; Athanas- 
soula et al. 2009a; Harsoula & Kalapotharakos 
2009; Athanassoula et al. 2009b, 2010; Contopou- 
los & Harsoula 2011). In this scenario the spiral 
dynamics are coupled to the bar, and are driven 



by the manifolds. This approach has been studied 
by two different groups of people. 

One of those groups (Romero-Gomez et al. 
2006, 2007; Athanassoula et al. 2009a,b, 2010), 
considers a continuous flow of orbits along the 
manifolds emanating from Li or L2. When spi- 
rals form, stars move away from the corotation 
in a radial movement (Athanassoula et al. 2010), 
and material is needed to replenish the mani- 
folds. One prediction of this "manifold theory" 
(or "Lyapunov tube model"), not accounted for in 
the density wave scenario, is that stronger bars 
should have more open spirals as compared to 
weaker bars, i.e., the spiral arms pitch angle^ 
should increase with bar strength (Athanassoula 
et al. 2009b). This kind of correlation was previ- 
ously predicted by Schwarz (1984), although for 
gas arms driven by a bar perturbation. 

Another view of the "invariant manifold the- 
ory" (Voglis et al. 2006b,c; Tsoutsis et al. 2008, 
2009) considers the locus of all points with initial 
conditions at the unstable manifolds that reach a 
local apocentric (or pericentric, see Harsoula et 



^The angle between a tangent to the spiral arm at a certain 
point and a circle, whose center coincides with the galaxy's, 
crossing the same point. 
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al. 2011) passage, i.e., the apsidal sections of the 
manifolds. In this scenario, there is no need for 
the replenishment of material to obtain long-lived 
spirals (see, e.g., Efthymiopoulos, C. 2010). Both 
views of the "invariant manifold theory" predict a 
trailing spiral pattern for strong perturbations and 
similar pattern speeds for the bar and spiral, i.e., 
r^bar = Qspirai^ Howcvcr, in the view of Voglis et 
al. (2006b,c) and Tsoutsis et al. (2008, 2009), the 
"azimuthal tilt" of the spiral response (Tsoutsis 
et al. 2009), i.e., the difference between the bar's 
major axis and the Lagrangian points Li or L2 at 
the moment of the onset of the spiral, determines 
how open the spiral arms will be. In this case, the 
pitch angles are smaller than the ones predicted 
by Athanassoula et al. (2009b) and become even 
smaller for pure bar models when the "azimuthal 
tilt" is not taken into account (C. Efthymiopoulos, 
private communication 2011). 

Patsis et al. (2010) describe one more dynami- 
cal mechanism that supports spiral arms through 
stars in chaotic motion. They propose this mech- 
anism by describing the spiral arms of the barred- 
spiral NGC 1300. Together with the bar, these 
spiral arms are inside the corotation and are not 
related to the presence of unstable Lagrangian 
points and the associated families of periodic or- 
bits. This alternative mechanism may be linked to 
some range of pitch angles of spiral arms encoun- 
tered in barred-spiral systems. 

Do manifolds drive spiral dynamics in barred 
galaxies? Or are the dynamics driven by the 
bar? The bar may drive the dynamics, affecting 
the spiral amplitude locally, as reported by Salo 
et al. (2010) (see also Block et al. 2004) and 
previously discarded (or weakly corroborated) by 
other authors comparing bar strength to spiral 
arm strength (Buta et al. 2009; Durbala et al. 
2009; Seigar & James 1998). Bars driving the 
dynamics would imply an accordance with (lin- 
ear) density wave theory. These spirals may be a 
continuation of the bar mode, or an independent 
mode coupled to the bar (e.g.. Tagger et al. 1987; 
Masset & Tagger 1997). In the "Lyapunov tube 
model", the strength of the bar affects the pitch 
angle of the spirals, but not its amplitude. The 
amplitude of the spirals depends on how much ma- 
terial is trapped by the manifolds, although, the 
amplitude of the spirals should in general decrease 
outward (Athanassoula et al. 2010). Grosb0l et 



al. (2004) investigated the relation between the 
amplitude of the spirals with the pitch angle in 
non-barred and weakly barred galaxies. 

One prediction of the density wave theory (Hozumi 
2003, see §6.1) entails that different pitch angles 
are expected for spirals when observed in different 
bands (e.g., optical versus near-infrared [NIR]). 
According to Athanassoula et al. (2010), the "in- 
variant manifold theory" predicts that stars of 
different ages will be guided by the same mani- 
fold, and no difference between the winding of the 
spirals is expected. 

In this paper, we will investigate whether the 
predictions of pitch angles are observed for real 
galaxies, or not. Two methods were applied for 
this purpose: the "slope method" (Section 4.1), 
which is especially good for determining how long 
the logarithmic shape is maintained for spiral 
arms, and the "Fourier method" (Section 4.2), 
which was used to determine the "dominant" pitch 
angle inside a given annulus for each object. 

2. GALAXY SAMPLE 



The initial galaxy sample consists of 104 galax- 
ies classified as Fourier bars in Laurikainen et 
al. (2004). The data were acquired from the 
Ohio State University Bright Galaxy Survey (OS- 
UBGS) (Eskridge et al. 2002). From this initial 
sample, it was found that only 84 objects present 
spiral-like features. Nevertheless, not all the ob- 
jects are suitable for this kind of study due to 
asymmetries, e.g., short, faint, or ragged spiral 
arms, or prominent rings. The following criteria 
were established in order to obtain a sample, in- 
cluding objects with a morphology candidate to 
be explained by "chaotic" spirals. 

1. The spiral arms must remain logarithmic, 
i.e., with a constant pitch angle (i), at least 
for 50° in the azimuthal range, a.^ This was 
verified with the "slope method" (see Sec- 
tion 4.1). The lower limit value of a was cho- 
sen according to Figure 4 in Athanassoula et 
al. (2009b), where the manifold loci remain 



^Although the spiral arms may extend further in the disk 
with a varying pitch angle, i.e., different slopes in a Inr 
versus 6 map. 
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logarithmic (for the adopted model param- 
eters) and maintain a "nearly" logarithmic 
geometry up to ~ 100°. We consider that 
the manifold loci and the density maximum 
along the spirals coincide. According to Pat- 
sis (2006), spirals supported by chaotic par- 
ticles may extend up to 7r/2 radians. Vari- 
ations of a toward larger angles will be dis- 
cussed in Section 6. 

2. The object presents two spiral arms visually 
connected to the bar. 

3. No prominent inner rings (near the bar's 
end) are present.^ Ring structures are con- 
nected to the bar on both sides. The pitch 
angle definition as applied in this investiga- 
tion only refers to spiral arms. A dependence 
of the inner ring shape on bar strength has 
been investigated by Grouchy et al. (2010). 

After applying these selection criteria, the final 
sample consists of 27 barred spirals (see, e.g.. Ta- 
ble 1). 

In order to use the bar strength values of Lau- 
rikainen et al. (2004, see Section 3), we adopt the 
same deprojection parameters of those authors, 
i.e. the same values for position angle (0) and 
minor-to-major axial ratios {q = b/a).^ To deter- 
mine these parameters, Laurikainen et al. (2004) 
fit ellipses to the outer isophotes on the disk. They 
were based on the OSUBGS 5-band images that 
are deeper than iJ-band images. 

To test the Athanassoula et al. (2009b, 2010) 
predictions regarding spiral arms pitch angles, we 
use the NIR i7-band since we are interested in 
"long" -lived structures rather than young stars, 
HIT regions, or gas that would be present in optical 
data. 



3. BAR STRENGTH 



^An exception is NGC 5921 where a ring is present, but it 
does not dominate over the spiral features. 

^With the exception of NGC 1300, for which we adopt (f) = 
100° ± 14, and q = 0.6± 0.1 (see Section 5). Uncertainties 
of 10% in incHnation translate in 10%-15% in perturbation 
strength (Buta & Block 2001; Laurikainen & Salo 2002). 



The predicted trend in Athanassoula et al.'s 
(2009a) "manifold models" requires the strength 
of the bar at the radius of the Lagrangian points 
Li or L2. It should be mentioned that for these 
models the self-gravity of the spirals was not taken 
into account. On the other hand, the addition of 
the spiral potential in Tsoutsis et al.'s (2009) mod- 
els shifts the positions of the Lagrangian points Li 
or L2 both in the radial and azimuthal directions. 

The strength of the bar can be obtained from 
the Laurikainen et al. (2004) radial profiles of the 
perturbation strength. Laurikainen et al. (2004) 
used the gravitational torque method (Combes & 
Sanders 1981; Buta & Block 2001; Block et al. 
2002) taking care of the artificial bulge stretch (see 
also Speltincx et al. 2008). The perturbation 
strength is calculated as 

Qtir) = (1) 

' dr 

which represents the ratio between the maximum 
amplitude (over azimuth) of the tangential force, 
and the mean axisymmetric radial force derived 
from the m = component of the gravitational 
potential. The potential is inferred from the lumi- 
nous mass, and can be represented as (see Combes 
& Sanders 1981; Quillen et al. 1994): 

^(r,6>) ^ ^o(^) + ^ ^rn{r) cos [mO]. (2) 

m=2,4,6 

The angle is given in the deprojected image, and 
^ = along the bar major axis. For this investi- 
gation we assume that Li = L2 = L. For real 
galaxies, Li may diff"er from L2 due to odd terms 
in the gravitational potential. 

We analyzed three cases in which the bar's 
strength is estimated in three different ways. 

1. In the first case, the bar's strength is esti- 
mated at r = tl- The Lagrangian point 
or corotation radius (Sellwood & Wilkinson 
1993), vl, was obtained from Buta & Zhang 
(2009), who applied the "potential- density 
phase shift method" to the OSUBGS sam- 
ple. There have been significant discussions 
on the validity of this method. This is partly 
because Zhang & Buta (2007) found some 
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cases (e.g., NGC 4665) where r^/rbar < 
1, i.e., corotation before the end of the 
bar. According to Contopoulos (1980), self- 
consistent bars are not possible to be mod- 
eled in this regime. 

One important difference between "mani- 
fold" models (Romero-Gomez et al. 2006; 
Voglis et al. 2006b, c; Athanassoula et al. 
2009a) and the "potential-density phase 
shift method" is that Zhang & Buta (2007) 
and Buta & Zhang (2009) considered po- 
tentials varying considerably with time. 
The time-independent (rigid) potentials of 
the "manifold" models generate "passive" 
chaotic orbit responses. Although the Zhang 
& Buta (2007) models involve chaos in the 
individual stars' trajectories, "collective dis- 
sipation" makes possible the existence of 
coherent structures (e.g., spiral arms). We 
define "rezog" as the corotation radius ob- 
tained from Buta & Zhang (2009). Table 1 
shows the Qt{r = tbzoo) values for the 27 
OSUBGS barred galaxies. 

2. In the second case, we estimate the bar's 
strength at a distance vl = 1.2rbar- Accord- 
ing to various studies (Athanassoula 1992; 
Elmegreen 1996; Aguerri et al. 2003), the ex- 
pected range for the bar length lies between 
^bar = tl/I-O and rbar = Elmegreen 
(1996) and Aguerri et al. (2003) also discuss 
objects where rbar = vl/U. A mean value 
of ^bar = is cxpcctcd for large sam- 
ples of galaxies. For bar strengths, the effect 
of having vl = l.Orbar or vl = 1.4rbar, in- 
stead of = 1.2rbar5 could be much larger 
than deprojecting a galaxy within 10% er- 
ror in the projection angles. For this study 
the bar length, rbar, was taken from Lau- 
rikainen et al. (2004). In Table 1 we show 
the Qt(r = 1.2rbar) values adopted for this 
investigation. 

3. The third case involves the maximum of the 
radial Qt{r) profiles or Qg. These were tab- 
ulated in Laurikainen et al. (2004). 

The adopted Qt{r) values from Laurikainen et 
al. (2004) were computed assuming a constant 
M/L ratio throughout the disk and an empirical 
correlation for the vertical scale- height (hz). Also, 



it is assumed that dark matter has little impact 
on the bar strength. For the Qt(r) error calcula- 
tion shown in Table 1 , the Qg error of Laurikainen 
et al. (2004) was summed in quadrature with the 
error inherent to digitization^ of the Qt(r) plots 
and the r^ (Buta & Zhang 2009) errors for the 
Qt{r = rezog) values. 

A technique for separating the gravitational 
torques of bars and spirals was developed by Buta 
et al. (2003, 2005). This technique separates the 
bar-hdisk image to obtain the bar strength (at 
the respective maximum of Qt(r)) unaffected by 
the spiral gravitational infiuence. Nevertheless, 
for the majority of barred galaxies in the OSUBGS 
sample, the bar strength, Q^, dominates over the 
spiral arm strength Qg (Durbala et al. 2009). Also, 
the correction of the spiral arms does not affect the 
tendencies for Qg in the Hubble sequence (Lau- 
rikainen et al. 2007). In either case, for this inves- 
tigation it is assumed that Qg ^ Qh^ and that the 
Qt{T) values are affected by the spirals within the 
errors. 

4. PITCH ANGLES 



Spiral arms pitch angles have been measured 
in the literature with different methods. Dan- 
ver (1942) measured the spiral arms on photo- 
graphic plates. Kennicutt (1981) measured the spi- 
ral shapes using the intensity and HII region dis- 
tributions. Ma et al. (1999) fit the shapes of spiral 
arms directly on the images. Fourier decomposi- 
tion methods had also been used (e.g., Considere 
& Athanassoula 1988; Puerari & Dottori 1992; 
Saraiva Schroeder et al. 1994; Seigar et al. 2006), 
yielding similar results as other methods (Con- 
sidere & Athanassoula 1988; Puerari & Dottori 
1992). 

4.1. "Slope Method" 

This method is similar to the one used in Seigar 
& James (1998). It is assumed that the arms can 
be represented by logarithmic spirals, which im- 
plies a constant pitch angle. Although, variable 
pitch angles may be a better and more adequate 



^The ADS's data extraction applet, DEXTER (Demleitner 
et al. 2001), was used for this purpose. 
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Fig. 1.— Deprojected i7-band images of NGC 150, NGC 210, NGC 289, NGC 578, NGC 864, NGC 1073, 
NGC 1187, and NGC 1300. The display is in a logarithmic scale. The analyzed arm segments for the "slope 
method" are shown in the figures (solid lines), together with the annulus adopted for the "Fourier method" 
(dashed lines). 
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Fig. 2.— Deprojected i7-band images of NGC 1703, NGC 1832, NGC 3059, NGC 3261, NGC 3513, NGC 
3583, NGC 3686, and NGC 4145. The display is in a logarithmic scale. 
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Fig. 3.— Deprojected iJ-band images of NGC 4303, NGC 4902, NGC 4930, NGC 4995, NGC 5483, NGC 
5921, NGC 6221, and NGC 6300. The display is in a logarithmic scale. 
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Fig. 4.— Deprojected i7-band images of NGC 6384, NGC 7479, and IC 5325. The display is in a logarithmic 
scale. 



Table 1 
Perturbation Strengths 



Galaxy 


Qt{r 


— rBzoo) 


Qt{r 


= 1.2rbar) 


Gah 


ixy 


Qt{r 


— ^Bzog) 


Qt{r 


= 1.2rbar) 


NGC 150 


0.302 


±28.1% 


0.352 


±24.1% 


NGC 


3686 


0.075 


±24.0% 


0.152 


±11.8% 


NGC 210 


0.060 


±5.0% 


0.058 


±3.4% 


NGC 


4145 


0.119 


±2.5% 


0.122 


±2.5% 


NGC 289 


0.133 


±3.8% 


0.099 


±7.1% 


NGC 


4303 


0.149 


±29.5% 


0.251 


±17.5% 


NGC 578 


0.092 


±79.3% 


0.040 


±37.5% 


NGC 


4902 


0.068 


±38.2% 


0.128 


±20.3% 


NGC 864 


0.201 


±18.9% 


0.236 


±15.7% 


NGC 


4930 


0.138 


±15.9% 


0.099 


±22.2% 


NGC 1073 


0.498 


±6.4% 


0.386 


±3.4% 


NGC 


4995 


0.280 


±16.8% 


0.263 


±17.9% 


NGC 1187 


0.162 


±26.5% 


0.206 


±20.9% 


NGC 


5483 


0.098 


±4.1% 


0.111 


±3.6% 


NGC 1300 


0.475 


±3.4% 


0.277 


±4.0% 


NGC 


5921 


0.384 


±6.0% 


0.329 


±7.0% 


NGC 1703 


0.100 


±5.0% 


0.058 


±8.6% 


NGC 


6221 


0.309 


±36.2% 


0.237 


±47.3% 


NGC 1832 


0.121 


±19.8% 


0.149 


±16.1% 


NGC 


6300 


0.158 


±6.3% 


0.080 


±3.8% 


NGC 3059 


0.251 


±20.3% 


0.268 


±17.9% 


NGC 


6384 


0.136 


±14.7% 


0.048 


±41.7% 


NGC 3261 


0.133 


±6.8% 


0.117 


±7.7% 


NGC 


7479 


0.516 


±11.6% 


0.240 


±25.0% 


NGC 3513 


0.256 


±27.0% 


0.210 


±32.9% 


IC 5325 


0.211 


±9.5% 


0.102 


±19.6% 


NGC 3583 


0.160 


±5.0% 


0.207 


±2.9% 















Note. — Columns (1) and (4): object name. Columns (2) and (5): perturbation strength (see Section 3) 
from Laurikainen et al. (2004), at corotation radius from Buta Sz Zhang (2009). Columns (3) and (6): 
perturbation strength from Laurikainen et al. (2004), at radius r = 1.2rbar (see Section 3). 
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representation for some objects (see e.g., the case 
of NGC 1365 in Ringermacher & Mead 2009). 

Before deprojection, the spiral regions were iso- 
lated by masking the bar, foreground stars, strong 
star- forming regions (visually selected), bad pix- 
els, and other structures not associated with the 
corresponding arm region. After deprojection {H- 
band data), the centers of the objects were deter- 
mined by fitting ellipses to the central isophotes 
close to the bar region.^ Afterward, the spiral 
arms were "unwrapped" by plotting them in a In r 
versus map (e.g., lye et al. 1982; Elmegreen et 
al. 1992; Grosb0l et al. 2004). Under this geomet- 
ric transformation, logarithmic spirals appear as 
straight lines. The pitch angle, i, is related to the 
slope of the line, s, as 

cot i = k\s\, (3) 

wherein A: is a constant^ due to "pixelation" and 
unit conversion. 

Two arm segments were selected closest to the 
bar's end with the condition that the slope, 5, was 
maintained nearly constant along them (see Fig- 
ures 1-4) . Due to this "slope restriction" , in many 
cases the critical segment including the part of the 
arms attached to the bar, was not able to be con- 
sidered. The slope (s) is determined by first se- 
lecting for each column in the arm segment the 
pixels with a maximum in intensity (see as an ex- 
ample Figure 5, for the case of NGC 1832). A 
least-squares fit is then obtained for the result- 
ing pixels. As already mentioned, these fits were 
done in the H-hand aiming to trace Population 
II stars. Young stars and clusters can contribute 
locally up to 20%-30% of the observed radiation 
in the NIR (e.g., Rix & Rieke 1993; Rhoads 1998; 
Patsis et al. 2001; Grosb0l & Dottori 2008). How 
these young objects affect the pitch angles' mea- 
surements depends on the star formation condi- 
tions and young stars kinematics. For this investi- 
gation, it is assumed that young stars and clusters 



^This was done with the ELLIPSE task in IRAF. IRAF is 
distributed by the National Optical Astronomy Observa- 
tory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agree- 
ment with the National Science Foundation. 

^Basically converts ^-axis pixels to radians, and determines 
the equivalence between pixels in the In r axis and physical 
units of an image. 



affect the spiral arms pitch angles within the errors 
involved in the methods applied. 

As previously mentioned, all the objects with 
inner rings, asymmetries, unclear, or "logarithmi- 
cally short" arms were discarded from the analysis. 
For the remaining 27 objects, the arm segments (I 
or II) best determined and with the clearest spiral 
structure were also identified. These are marked 
with an asterisk (*) in Table 2, together with the 
adopted radial ranges, Ar, tabulated from inner- 
most (ro) to outermost radius (r), and azimuthal 
ranges, a. Azimuthal ranges are obtained by the 
equation 

a = cot ii^ In ( — ) (4) 

and are displayed graphically in Figures 1-4 with 
regions delimited by solid lines. These values do 
not indicate the end of the spirals, since spiral 
arms may extend further with a variable pitch 
angle (Ringermacher & Mead 2009). Although if 
the extensions of the spirals have reduced ampli- 
tudes with respect to the logarithmic part, their 
phases will be difficult to determine, as will their 
pitch angles. The estimations are done indepen- 
dently of the amplitude (strength) of the spiral 
itself. Grosb0l et al. (2004, their Figure 8) found 
a tendency between the amplitude of the m = 2 
spiral and pitch angles in SA and SAB galaxies. 

Figure 6 shows a histogram of the maximum 
azimuthal range distribution (either arm segment 
I or II) for each object presented in Table 2. 

4.1.1. Error Determination 

Errors introduced by deprojection parameters 
(0 and q) translate into different slopes or devia- 
tions of a straight line in the Inr versus plots. 
For each object, five deprojected frames were ob- 
tained to better account for these errors. The im- 
ages were deprojected with the parameters ^, g; 

sd^q; (j) — sd, q; (j)^q^ sd; and 0, g — sd^ where 
sd is the respective standard deviation. Pitch an- 
gle values were measured and compared to the case 
when (j) and q were used as the deprojection pa- 
rameters (i.e., when sd = 0). The cases with the 
highest (positive) or lowest (negative) discrepan- 
cies were adopted to account for the +(7 and —a 
errors, respectively (see Table 2). 
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Table 2 

Slope Method" Derived Parameters 



Galaxy and Segment if/ (deg) Ar(arcsec) Q;(deg) Segment if/ (deg) Ar(arcsec) Q;(deg) 



NGC 150 Arm 


24 


£3 + 5.3 
6-1.5 


(34.7-55.4) 


58 ± 8 


Arm II 


qq 9+7.6 


(41.8-51.4) 


18 ± 5 


NGC 210 Arm I* 


18 


A+0.7 
^-0.9 


(64.8-107.6) 


88 ± 4 


Arm II 




(63.6-105.6) 


51 ± 4 


NGC 289 Arm I 


16 


^-1.2 


(22.8-28.2) 


41 lb 4 


Arm II* 


^ „ 0+4.8 

19-3Io.9 


(22.8-30.8) 


50 ± 7 


NGC 578 Arm I* 


23 


Q + 0.4 
^-1.3 


(20.8-50.3) 


115 ± 4 


Arm II 




(20.8-50.3) 


110 ± 7 


NGC 864 Arm I* 


24 


q + 0.7 


(33.0-68.5) 


93 ± 3 


Arm II 


OO 0+4.1 

28.81^ 2 


(34.8-53.1) 


44 ± 4 


NGC 1073 Arm 1^ 


29 


n + 7.8 
0-6.7 


(37.8-55.3) 


39 ± 10 


Arm II 


12.6^1 8 


(48.7-70.0) 


93 ± 18 


NGC 1187 Arm I 


35 


7.+3.8 
' -4.1 


(35.7-48.8) 


25 ± 3 


Arm II* 


29-lto.5 


(25.2-55.6) 


82 ± 6 


NGC 1300 Arm I* 


21 


^-7.8 


(84.8-116.9) 


47 db 14 


Arm II 


18.8t^;^ 


(74.3-100.5) 


51 ± 13 


NGC 1703 Arm I* 


16 


7+0.4 
— 1.9 


(15.6-28.1) 


112 ± 7 


Arm II 


25.ot°:7 


(16.8-28.1) 


64 ±2 


NGC 1832 Arm I* 


20 


7+1.3 
' -1.2 


(20.0-45.8) 


126 ±8 


Arm II 


27.01^-9 


(20.0-27.8) 


37 ±6 


NGC 3059 Arm I* 


31 


c: + l-9 
^-1.6 


(26.8-54.0) 


66 ±4 


Arm II 


33.7l^;5 


(23.6-33.9) 


31 ±5 


NGC 3261 Arm I* 


17 


-,+15.6 
-^-2.5 


(27.4-34.2) 


41 ± 14 


Arm II 




(30.0-36.1) 


50 ± 7 


NGC 3513 Arm I 


32 


Q+1.3 

^-0.9 


(26.1-72.6) 


91 ±4 


Arm II* 


40.7ti;^ 


(16.3-72.6) 


100 ±5 


NGC 3583 Arm I 


20 


4+0.5 
^-3.3 


(34.7-59.3) 


83 ± 7 


Arm II* 


24 1+2-° 


(34.7-68.8) 


87 ±5 


NGC 3686 Arm I* 


24 


g + 0.9 
"-0.5 


(20.2-46.2) 


104 ±3 


Arm II 




(20.2-43.7) 


19 ±4 


NGC 4145 Arm I* 


17 


^-0.4 


(50.5-87.9) 


104 ±8 


Arm II 




(50.5-87.9) 


100 ± 12 


NGC 4303 Arm I* 


38 


4+10 

^-1.2 


(29.2-60.5) 


53 ±2 


Arm II 


44.111- 


(29.2-57.3) 


40 ±2 


NGC 4902 Arm I* 


11 


9 + 6.6 
^-1.8 


(21.4-30.0) 


98 ±26 


Arm II 


OQ 1 +5.3 


(19.5-34.7) 


77 ±9 


NGC 4930 Arm I* 


18 


7+0.6 
' -2.7 


(38.9-58.2) 


68 ±6 


Arm II 


74 A + 5-5 
'^•"-16.1 


(36.2-53.1) 


6±4 


NGC 4995 Arm I 


47 


+ 3.7 
^-4.9 


(22.6-34.3) 


22 ± 3 


Arm II* 




(17.3-37.7) 


56 ±8 


NGC 5483 Arm I* 


29 


9 + 0.5 
^-4.4 


(14.8-39.4) 


100 ± 9 


Arm II 


31.813^ 


(14.8-34.7) 


79 ± 6 


NGC 5921 Arm I* 


16 


1+0.5 
-^-0.6 


(61.0-95.3) 


89 ± 3 


Arm II 




(69.5-91.9) 


54 ± 3 


NGC 6221 Arm I* 


19 


r + 1-3 
^-0.7 


(39.0-57.6) 


63 ± 3 


Arm II 




(41.2-56.5) 


50 ± 4 


NGC 6300 Arm I* 


17 


O + 0.6 
"^-0.4 


(42.2-78.9) 


116 ± 3 


Arm II 


i8.3tS:t 


(42.2-78.9) 


108 ± 2 


NGC 6384 Arm I* 


24 


cc + 0.4 
^-1.9 


(31.4-84.9) 


125 ± 6 


Arm II 




(40.1-57.2) 


63 ± 4 


NGC 7479 Arm I* 


36 


g + 2.3 
"-1.4 


(56.1-95.8) 


41 ± 3 


Arm II 


23.0to^ 


(60.4-95.8) 


62 ± 5 


IC 5325 Arm I 


49 


9 + 2.3 
^-9.7 


(11.4-23.0) 


35 ± 7 


Arm II* 




(14.3-27.4) 


100 ± 8 



Note. — Column (1): object and spiral arm segment, see Figures 1-4. Columns (2) and (6): ii"-band pitch angles, 
in, in degrees. Columns (3) and (7): radial ranges, ro to r, in arcsec. Columns (4) and (8): azimuthal ranges , 

a = cot ijf In f ^ j , in degrees. Column (5): spiral arm segment (see Figures 1-4) for the same object as Column (1). 
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4.2. "Fourier Method" 



Figure 7(a) plots the pitch angles in arm seg- 
ment I versus arm segment II for each object as 
obtained with the "slope method". Figure 7(b) 
shows a histogram of the absolute value differ- 
ence between arm segments I and 11. As shown in 
the figures, some scatter is present when analyz- 
ing spiral arm segments within the same galaxies. 
Since we are interested in comparing single values 
of pitch angles for each object, we need a method 
that provides the "dominant mode" for the pitch 
angle measurement. The "Fourier method" is per- 
fectly adequate for this purpose. 

In this method, it is again assumed that the 
arms can be represented by logarithmic spirals.^ 
The Fourier amplitudes for each component are 
given by 



A{m,p) 



J2i=i Ej=i r, 0)exp[-i{mO + p\m 
ELiE/=i/^.(lnr,^) 

(5) 

where r and 6 are the polar coordinates, lij is the 
intensity at coordinates Inr, 0, m is the number 
of spiral arms (or modes), and p is related to the 
spiral arms pitch angle (in) by 



-m/pn 



(6) 



where Pmax corresponds to the maximum of 
A{m,p) and m = 0, 1, 2, 3, . . ., i.e., the maximum 
of the Fourier spectrum (see, e.g., Puerari & 
Dottori 1992; Saraiva Schroeder et al. 1994) for 
mode m. Most of the analyzed objects present 
m = 2 as the dominant mode for the spiral arms 
in the i7-band (see Table 3), so it was adopted for 
this investigation. The exceptions are NGC 3261 
and NGC 4930 in which m = 1 dominates and was 
used instead. For NGC 1300 and NGC 7479, other 
Fourier modes (m) compete with the m = 2 mode 
because of the spiral arm segments with variable 
pitch angles. The pitch angles corresponding to 
the m = 2 Fourier mode were adopted for these 
objects in the subsequent analysis. 

For the galaxies of the sample, it has been real- 
ized that the presence of foreground stars does not 



^However, a Fourier analysis can be done without the as- 
sumption of a constant pitch angle. 



affect the value of the pitch angle in general. Nev- 
ertheless, caution must be taken when foreground 
stars (or objects) compete in extension with spiral 
arms (see, e.g., annulus for NGC 864 in Figure 1). 
In these cases the need for masks is required. 

Objects were deprojected as explained in Sec- 
tion 4.1. Radial ranges were selected to cover 
the spiral segments previously analyzed with the 
"slope method" . The azimuthal coverage is 27r ra- 
dians. The analyzed annuli are shown graphically 
in Figures 1-4 (dashed lines). These are the re- 
gions where the Fourier analysis was performed. 

Table 3 shows the results for the Fourier pitch 
angle values, which agree with the "slope method" 
within a ^ 16° difference (this corresponds to lo- 
in Figure 7) in the majority of the objects. NGC 
5921 and NGC 6221 present the largest differences 
(- 15°). For two objects, NGC 4995 and IC 5325, 
the computed pitch angles are close to ~ 90° . This 
is due to the fact that the spiral arms have a low 
surface brightness (as compared to the disk) and 
the bar component is difficult to isolate in the ana- 
lyzed annulus. The "slope method", for the "best- 
defined arm", was used instead for these two ob- 
jects in the subsequent analysis. 

4-2.1. Error Determination 

Errors were determined in the same way as in 
the "slope method" . These were added in quadra- 
ture with the error intrinsic to the method. A pro- 
gram was built that computes the two-dimensional 
fast Fourier transform in Equation 5. The output 
of this program is a 128 x 2048 {m,p) matrix. The 
two closest values near Pmax were used to approx- 
imate the error of the method. 



5. COMMENTS FOR SOME OBJECTS 



NGC 210. "Skinny" spiral arms compared with 
the bar. 

NGC 289. The outer spiral arms have a greater 
pitch angle {in ~ 40°) as compared to the inner 
ones {in ^ 25°; "slope method"). 

NGC 578. Two symmetric spiral arms near the 
bar. 
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Azimuthal range (degrees) 



Fig. 6. — Histogram of the maximum azimuthal ranges (either arm segment I or II, see Table 2). 




(b) 



n n n 



20 40 60 

Pitch angle arm I (deg) 



10 20 30 40 50 60 
Difference between arm I & arm II (degrees) 



Fig. 7. — (a) Pitch angles (in deg) for arm segments I (x-axis) vs. arm segments II (?/-axis) for each object. 
Dotted line: one-to-one relation, (b) Histogram of the absolute value difference between arm segments I and 
II, obtained for each object with the "slope method". The standard deviation around the zero difference 
value is 16? 3. 



12 




ln(r) 



Fig. 5. — Plot of Inr vs. for arm segment I in 
NGC 1832 (i^-band). Crosses indicate the points 
where a maximum intensity was found for each 
column for the corresponding section in the "un- 
wrapped" image. The continuous line indicates 
the least-squares fit. 



NGC 1073. Spiral arms difficult to trace (low 
signal-to- noise ratio). 

NGC 1187. The two arm features analyzed are 
visually attached to the bar. A third arm feature, 
not visually attached to the bar, is present. The 
radial ranges for the "Fourier method" were mod- 
ified with respect to the "slope method" to allow 
a better signal-to-noise ratio in the Inr versus 
map. 

NGC 1300. Two well-defined logarithmic spi- 
ral arms, although short in azimuthal range. The 
adopted deprojection parameters were changed as 
compared to the ones of Laurikainen et al. (2004). 
This was done because the values provided in Lau- 
rikainen et al. (2004) do not agree with the outer 
isophotes of the OSUBGS images. An average be- 
tween Hyperleda (Paturel et al. 2003), RC3 (de 
Vaucouleurs et al. 1991), and a visual determina- 
tion of the outer isophotes was used. 

The deprojection parameters from Lindblad et 
al. (1997) were also tried for the pitch angle mea- 
surements. These parameters, (j) = 87° ± 2° and 
q = 0.82 ±0.05, are based on H I data and are inde- 
pendent of kinematical or dynamical criteria (see 
also Kalapotharakos et al. 2010). Using these pa- 
rameters, spiral arms are difficult to follow in a 
Inr versus map (assuming a logarithmic geome- 
try). For arm region I, a pitch angle of 21? 2 ± 6° 
was obtained. Arm region II was not possible to 
measure via the "slope" method. The pitch angles 
obtained by applying the "Fourier" method led to 
values with a contrary sign to the one expected, 
i.e., an inverse sense of winding for the spiral arms. 

NGC 1703. Difficult to analyze the spiral arms 
in the inner regions due to few pixels in a Inr 
versus map. 

NGC 1832. The bar region is distorted (not 
straight). 

NGC 3059. "Hard to follow" logarithmic shape 
for the spiral arms. 

NGC 3583. Two symmetric spiral arms can be 
appreciated in the outer disk. The region close to 
the bar presents a structure similar to a ring or a 
tight spiral arm. 

NGC 414^- Double bar system? 

NGC 4303. This object presents three main 
spiral arms. 

NGC 4902. Three spiral regions are present in 
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this object. 

NGC 5921. This object presents an inner ring 
and spiral features. 

NGC 6300. This object presents spiral features 
and apparently a ring feature. 

NGC 6384' Spiral arms with bifurcations. 

NGC 7479. In general the spiral arms for this 
object do not present a clear logarithmic geometry. 

IC 5325. This object presents four well defined 
segments of spiral arms. Only the ones near the 
bar's end were analyzed. 

6. RESULTS AND DISCUSSION 



Figure 8 shows the results for the pitch an- 
gle, iH (Fourier method, except for NGC 4995, 
and IC 5325, see Section 4.2), versus perturba- 
tion strengths, Qtir = rezog)- A first inspec- 
tion of the data, where the "azimuthal range" ^ 
is a > 50°, shows considerable scatter around the 
predicted correlation for models A (Ferrers 1877, 
bar potential) and D (Dehnen 2000, bar potential) 
in Athanassoula et al. (2009b). However, if the 
a criterion is changed to logarithmic spiral seg- 
ments that extend up to a > 70°, a > 90°, and 
a > 110°, the scatter is reduced. The reduced 
Pearson's chi-square, X^/^^ obtained as 



where ik is the kth Fourier- measured pitch angle 
and ip is the predicted pitch angle value for models 
A and D in Athanassoula et al. (2009b), gives the 
results 3.10, 1.55, 1.83, and 2.00 for a > 50° (n = 
27), a > 70° (n = 17), a > 90° (n = 13), and 
a > 110° (n = 5), respectively. 

Figures 9 and 10 show the results for the cases 
Qt{T = 1.2rbar) and Q^, respectively. For a > 70°, 



^This is obtained with the "slope method" via Equation 4. 
It is the "maximum" azimuthal range that is taken into 
account, i.e., the greatest value of a for either arm segment 
I or 11. 

third model with a Barbanis & Woltjer (1967) bar po- 
tential (BW model) was considered in Athanassoula et 
al. (2009a,b). This model agrees with model D up to 
Qt(r = tl) ~ 0.2, and deviates toward higher pitch an- 
gles afterward, up to ~ 5° Qt{r = tl) ^ 0.6. 



reduced Pearson's chi-square values obtained as 

^2^j^ {Qk-Q,? ^ (8) 
k=i 

where Q/c is the kih. bar strength value correspond- 
ing to the kih Fourier-measured pitch angle and 
Qp is the predicted bar strength value for models A 
and D in Athanassoula et al. (2009b), yield the re- 
sults 0.049, 0.075, and 0.084 for the Qt{r = tezoq), 
Qt{r = 1.2rbar), and Qg plots, respectively. 

According to this result, the best concordance 
with the Athanassoula et al. (2009b) model is ob- 
tained by comparing the pitch angles with Qti^) 
given at Buta & Zhang (2009) bar corotation radii 
(r = rBzog)- This last point is not discussed 
in Athanassoula et al. (2009b). 

One important aspect in the Athanassoula et 
al. (2009b) prediction is that the self- gravity of 
the spirals was not taken into account. The poten- 
tial created by the "confined" chaotic orbits is ne- 
glected. Contrarily, Tsoutsis et al. (2009) empha- 
size the contribution of the spiral part for studying 
the dynamics of the "chaotic" spirals. Also, realis- 
tic bar potentials are hard to model. If many dif- 
ferent realistic potentials are used, the predicted 
correlations may become broader (Athanassoula 
et al. 2009b, 2010). This may explain in Fig- 
ure 8 the tendency of the points (squares and 
circles) to be above the predicted correlation for 
Qt(r = rBzo9) <0.2. 

6.1. Density Wave Theory Prediction 

The modal approach explains the density wave 
phenomena as generated by intrinsic mechanisms 
in the disk (Bertin et al. 1989a, b; Bertin & Lin 
1996). Normal modes of oscillation generate spon- 
taneously and evolve according to the physical and 
dynamic properties of the system. Three physi- 
cal properties determine the morphology in disk 
galaxies: the disk mass, the gas content, and the 
stellar velocity dispersion. When the disk mass is 
"high" , bar structures are generated as oscillating 
modes of the system. The modal theory considers 
bars and spirals equally, i.e., as normal modes of 
oscillation in the disk. 

Based on the dispersion relation, linear density 
wave theory predicts (Hozumi 2003) that the pitch 



14 



Table 3 

Fourier Method" Derived Parameters 



Galaxy z/f (deg) is (deg) Ar (arcsec) 



NGC 


150 





.416 


1 


.000 


0. 


,344 


0, 


.213 


0. 


,395 


0. 


,167 


27, 


Q + 0.9 
— 1.9 




(34.7-55.4) 


NGC 


210 





.525 


1 


.000 


0. 


,189 


0, 


.478 


0. 


,255 


0. 


,220 


16, 


7+0.5 
' —0.5 




(63.6-107.6) 


NGC 


289 





.103 


1 


.000 


0. 


,309 


0, 


.213 


0. 


,091 


0. 


,109 


19, 


7+0.8 
' —0.5 


17.2i;-| 


(22.8-30.8) 


NGC 
NGC 


578 
864 






.656 
.721 


1 
1 


.000 
.000 


0. 
0. 


,289 
,462 


0. 
0, 


.255 
.361 


0. 
0. 


,080 
,307 


0. 
0. 


,116 
,325 


24, 
20, 


9 + 1.6 

— 1.5 
7+1.2 

— 1.5 


18.01°;^ 


(20.8-50.3) 
(33.0-68.5) 


NGC 


1073 


0, 


.778 


1 


.000 


0. 


,515 


0, 


.506 


0. 


,307 


0. 


,349 


34. 


7+1.4 
' —5.9 




(37.8-70.0) 


NGC 
NGC 


1187 
1300 


0, 
0, 


.668 
.617 


1 
1 


.000 
.000 


0, 

1. 


,460 
,001 


0. 
0, 


.436 
.854 


0, 
0. 


,655 
,399 


0, 
0. 


,239 
,330 


21, 
11, 


4+1-3 

— 1.2 
9 + 15.3 

— 0.3 


13.11^0:3 


(29.0-59.4) 
(74.3-116.9) 


NGC 


1703 


0, 


.245 


1 


.000 


0. 


,255 


0, 


.371 


0. 


,161 


0. 


,196 


18, 


+ 1.5 
■^-2.0 


17.2lJ;« 


(15.6-28.1) 


NGC 


1832 


0, 


.242 


1 


.000 


0, 


,492 


0. 


.399 


0, 


,100 


0, 


,128 


25, 


1+1.8 
-^-1.1 




(20.0-45.8) 


NGC 


3059 





.728 


1 


.000 


0. 


,371 


0, 


.475 


0. 


,424 


0. 


,694 


27, 


■^-2.6 




(23.6-54.0) 


NGC 


3261 


1 


.043 


1 


.000 


0. 


,246 


0, 


.226 


0. 


,259 


0. 


,322 


9.1 


+ 13.3 
-0.3 


10.2^1",' 


(27.4-36.1) 


NGC 


3513 


0, 


.481 


1 


.000 


0, 


,229 


0. 


.472 


0, 


,141 


0, 


,220 


25, 


7+1.2 
• ' -1.7 




(27.8-72.6) 


NGC 


3583 





.446 


1 


.000 


0. 


,470 


0, 


.225 


0. 


,300 


0. 


,146 


24, 


r: + 2.4 
■^-2.1 


22.6l^;0 


(34.7-68.8) 


NGC 


3686 


0, 


.450 


1 


.000 


0. 


,886 


0, 


.271 


0. 


,276 


0. 


,493 


14, 


4+0-4 

■^-0.6 


14.4lS-« 


(20.2-46.2) 


NGC 


4145 


0, 


.448 


1 


.000 


0. 


,250 


0, 


.348 


0. 


,240 


0. 


,275 


23, 


g+1.0 

■^-0.7 


17.3l°-i 


(50.5-87.9) 


NGC 


4303 





.260 


1 


.000 


0, 


,174 


0. 


,294 


0, 


,149 


0, 


,186 


42, 


■^-1.9 




(29.2-60.5) 


NGC 
NGC 


4902 
4930 


0, 

1, 


.467 
.081 


1 
1 


.000 
.000 


0. 
0. 


,210 
,550 


0, 
0, 


.665 
.272 


0. 
0. 


,277 
,181 


0. 
0. 


,401 
,223 


20, 
30. 


a + 2.5 
■^-1.6 
1+2.1 
-^-2.1 


25.3tll 

i3.9t°:? 


(19.5-34.7) 
(36.2-58.2) 


NGC 


4995 





.898 


1 


.000 


0, 


,356 


0. 


,595 


0, 


,234 


0, 


,204 


90, 


Q + 5.6 
•'-'-8.8 




(22.5-37.7) 


NGC 


5483 


0, 


.163 


1 


.000 


0. 


,134 


0, 


.439 


0. 


,158 


0. 


,280 


26, 


+ 1.9 
■"^-1.8 




(14.8-39.4) 


NGC 


5921 


0, 


.188 


1 


.000 


0. 


,454 


0, 


.524 


0. 


,152 


0. 


,313 


30, 


g+4.9 
■"-3.1 


22.5tl° 


(61.0-95.3) 


NGC 


6221 


0, 


.602 


1 


.000 


0, 


,315 


0. 


.413 


0, 


,115 


0, 


,235 


35, 


4 + 4.7 
■^-1.4 




(39.0-57.6) 


NGC 


6300 





.199 


1 


.000 


0. 


,201 


0, 


.330 


0. 


,157 


0. 


,283 


23, 


1 +0.7 
■-^-0.9 


20.315;^ 


(42.2-78.9) 


NGC 


6384 


0, 


.501 


1 


.000 


0. 


,574 


0, 


.312 


0. 


,385 


0. 


,273 


26, 


+ 1.2 
"^-1.2 


28.0+_ro 


(31.4-84.9) 


NGC 


7479 


1, 


.073 


1 


.000 


0. 


,656 


1, 


.200 


0. 


,610 


0. 


,473 


26, 


+ 3.8 
•^-1.8 


'^-'-•^-3.2 


(56.1-95.8) 


IC 5325 





.367 


1 


.000 


0. 


,426 


0, 


.176 


0. 


,200 


0. 


,090 


69. 


4+3.8 
■^-5.2 


Sl-S-12^7 


(11.4-27.4) 



Note. — Column (1): galaxy name. Column (2): ratio between the maximum amplitudes of Fourier modes 
m = 1 and m = 2, in the if-band. Columns (3), (4), (5), (6) and (7): ratio between the maximum 
amplitudes of the respective Fourier modes, in the //-band. Column (8): //-band (see Section 4.2) pitch 
angles, in degrees. Column (9): S-band (see Section 6.1) pitch angles, in degrees. Column (10): radial 
ranges, ro to r, in arcsec. 
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T3 



20 




• a > 50 degrees 

A a > 70 deg. 

□ a > 90 deg. 

O « > 110 deg. 



0.2 0.4 0.6 

Qt (i" = ^'bzo9) 

Fig. 8. — Pitch angle in in deg vs. perturba- 
tion strength Qt{r = rezog) foi" the 27 galaxies 
selected for analysis (see Section 2). Lagrangian 
radius, = ^Bzog, from Buta & Zhang (2009). 
The dashed line corresponds to the predicted cor- 
relation for models A and D in Athanassoula et 
al. (2009b). Data are separated by a > 50° (ah 
points), a > 70° (triangles), a > 90° (squares), 
and a > 110° (circles). 
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Fig. 9. — Same as Figure 8 for perturbation 
strength Qt(r = 1.2rbar). 
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Fig. 10. — Same as Figure 8 for perturbation 
strength Qg. 
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angle should increase with increasing velocity dis- 
persion, or that 

taniocc^/S, (9) 

where i is the arms pitch angle, is the radial 
velocity dispersion, and H is the surface density 
of the disk. Spiral structure shows different mor- 
phologies when observed in optical versus NIR 
bands (Block & Wainscoat 1991; Grosb0l & Patsis 
1998). NIR bands can trace both the old popula- 
tions of bar and spiral arms, assuming that red 
young stars do not contribute globally to the ob- 
served radiation (Rix & Rieke 1993). Also, older 
populations have a higher velocity dispersion com- 
pared to younger ones (Barbanis & Woltjer 1967; 
Wielen 1977; Nordstrom et al. 2004; Binney & 
Tremaine 2008, section 8.4). For most galaxies at 
the arm location, we have that ^ 98% (by mass) 
of the stars belong to evolved populations (see 
Gonzalez & Graham 1996, and references therein). 
Nevertheless, young stars contribute to most of the 
light in optical wavelengths. According to this, 
NIR images of spiral perturbations should present 
higher pitch angles compared to optical ones. Az- 
imuthal age (color) gradients (e.g., Gonzalez 
& Graham 1996; Martmez-Garcia et al. 2009a,b; 
Martmez-Garcia & Gonzalez-Lopezlira 2011) may 
also affect the pitch angles observed in the op- 
tical versus NIR bands, but these are very diffi- 
cult to trace by just comparing the light distri- 
butions in two bands (Gonzalez & Graham 1996; 
Seigar & James 1998). Besides, azimuthal gradi- 
ents are not located continuously along the spi- 
ral arms but in specific regions (Gonzalez & Gra- 
ham 1996; Martmez-Garcia et al. 2009a; Martmez- 
Garcia & Gonzalez-Lopezlira 2011). 

From Equation 9, taking into account that 
young and old stars are similarly affected by the 
gravitational potential of the disk (which depends 
on the surface density), we obtain 

is = arctan I ^^^^ tanii:f|>, (10) 

where is the 5-band pitch angle, in is the H- 
band pitch angle, c^^ is the radial velocity disper- 
sion of young stars, and Crjj is the radial velocity 
dispersion of old stars. 

In the case of the invariant manifold theory. 



where chaotic orbits are "confined" in the spiral 
locus, no difference between pitch angles of spi- 
ral arms traced in different wavelengths is pre- 
dicted (Athanassoula et al. 2010). 

Seigar et al. (2006) found a nearly 1:1 correla- 
tion between pitch angle measurements in the B 
and H bands, for 57 galaxies in the OSUBGS (Es- 
kridge et al. 2002) sample. Nevertheless, based on 
the sample of five non-barred and weakly barred 
spirals, Grosb0l & Patsis (1998) notice that the 
main two-armed spiral is tighter when measured 
in bluer colors. For the barred-spirals data pre- 
sented in this investigation, we measured the pitch 
angles in the B-hand images for the same ob- 
jects analyzed in the //"-band from the OSUBGS 
sample, applying the "Fourier" method. The B- 
band images were registered to the iJ-band im- 
ages, so the high-resolution data (B-band) were 
degraded to the low-resolution data (i^-band in 
this case). Annulus regions were selected in the 
same positions as the iJ-band, and the pitch an- 
gles were measured identically^^ with the method 
described in Section 4.2. The results are shown^^ 
in Figure 11 (and Table 3) where a tendency of 
^ 30% of the points toward higher iJ-band pitch 
angles is observed. Although, if we apply the same 
azimuthal range (a) criteria as in Figures 8, 9, 
and 10, we can notice that ~ 80% of the a > 70° 
data lie very close to the 1:1 relation as expected 
(independently of a) from Athanassoula et al. 
(2010). 

6.2. Invariant Manifolds as Apsidal Sec- 
tions 

A comparison of the different treatments of the 
manifolds viewed as apsidal sections (Voglis et al. 
2006b,c; Tsoutsis et al. 2008, 2009), or as tubes 
that guide chaotic orbits (Romero-Gomez et al. 
2006, 2007; Athanassoula et al. 2009a,b, 2010), re- 
quires a different analysis involving separating spi- 
rals and bars. This will be covered in a subsequent 
publication. 



The same dominant modes, m = 1 or m = 2 (see Sec- 
tion 4.2), as measured in the if-band were adopted for the 
i^-band pitch angle measurements. 

'NGC 4995 and IC 5325 were excluded from this analysis 
(see the last paragraph in Section 4.2). 
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7. CONCLUSIONS 



The results of this investigation show the fol- 
lowing. 
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Fig. 11. — i^-band pitch angle vs. 5-band pitch 
angle, obtained via the "Fourier" method. Dotted 
line: one-to-one relation; short-dashed line: den- 
sity wave theory prediction from Equation 10, as- 
suming Cr^ ~ \f2crj^\ long-dashed line: density 
wave theory prediction from Equation 10, assum- 
ing Cr^ - 2C^^ . 



1. Although the adopted deprojection param- 
eters may introduce some biases (see, e.g., 
Barnes & Sellwood 2003), a trend can be 
observed where some strong barred spirals 
have more open spiral arms when compared 
to galaxies with weaker bars. This kind 
of trend was also discussed in Block et al. 
(2004), where a similar behavior was found. 
The correlation predicted by the manifold 
models of Romero-Gomez et al. (2006, 2007) 
and Athanassoula et al. (2009a,b, 2010) is 
better reproduced by observations on two 
conditions. 

(a) The corotation values obtained with 
the "potential-density phase shift method" 
& Zhang 2009) are adopted. 

(b) The spirals logarithmic geometry is 
maintained for large azimuthal ranges, 
OL > 70°. 

2. The ~ 60% of the 27 galaxies on the ana- 
lyzed sample seem to reproduce the investi- 
gated correlation. 

3. The pitch angles calculated via the "Fourier 
method" in the B (young stars) and the H 
(mostly old stars) bands yield similar values 
for ~ 80% of the objects where the azimuthal 
range, a, is greater than 70°. This kind of 
behavior is expected in the "Lyapunov tube 
model" (Athanassoula et al. 2010), although 
no restriction on the azimuthal range was 
given by the authors. 

4. Other possible mechanisms to generate spi- 
ral features in barred galaxies, such as bar- 
driven spirals (e.g. Salo et al. 2010), models 
where the Lagrangian points of the system 
are specified by both bar and spirals (e.g., 
Tsoutsis et al. 2009), or chaotic spirals in- 
side corotation (thus not related with the 
presence of unstable Lagrangian points; Pat- 
sis et al. 2010), cannot be excluded by the 
present investigation. 
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